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Abstract 

This  paper  presents  a  two-part  controller  that  uses  a  state 
variable  estimator  for  control  of  a  single  degree  of  freedom 
rotary  manipulator  actuated  by  Shape  Memory  Alloy  (SMA) 
wire.  A  model  for  the  SMA  actuated  manipulator  is  presented. 
The  model  includes  nonlinear  dynamics  of  the  manipulator, 
a  constitutive  model  of  the  Shape  Memory  Alloy,  and  the 
electrical  and  heat  transfer  behavior  of  SMA  wire.  The  es¬ 
timator  predicts  the  state  vector  at  each  time  step  and  corrects 
its  prediction  based  on  the  angular  position  measurements. 
The  transformation  temperatures  of  the  SMA  wire  are  stress 
dependent:  this  adds  to  complexity  of  the  control  of  the 
position  of  the  arm  for  certain  angular  positions.  To  overcome 
this  problem  a  sliding  mode  controller  is  designed  to  calculate 
the  desired  stress  of  the  wire  based  on  the  desired  angular 
position  of  the  arm.  A  classical  feedback  controller  maintains 
this  desired  stress  by  regulating  the  input  voltage  to  the  SMA 
wire.  The  stress  of  the  wire  is  a  better  set-point  for  the 
controller  since  the  transformation  temperatures  and  hence 
the  Martensite  fraction  are  functions  of  the  wire’s  stress. 
The  stabilization  and  tracking  performance  of  this  two-part 
controller  is  presented  and  compared  with  a  PID  controller. 

I.  Introduction 

SMAs  consist  of  a  group  of  metallic  materials  that  demon¬ 
strate  the  ability  to  return  to  some  previously  defined  shape 
or  size  when  subjected  to  the  appropriate  thermal  procedure. 
The  shape  memory  effect  is  hysteretic  and  occurs  due  to  a 
temperature  and  stress  dependent  shift  in  the  materials  crys¬ 
talline  structure  between  two  different  phases  called  Marten¬ 
site  and  Austenite  which  are  the  low  and  high  temperature 
phases,  respectively.  SMA  actuators  have  several  advantages 
for  miniaturization  such  as  excellent  power  to  mass  ratio, 
maintainability,  reliability,  and  clean  and  silent  actuation.  The 
disadvantages  are  low  energy  efficiency  due  to  conversion 
of  heat  to  mechanical  energy  and  inaccurate  motion  control 
due  to  hysteresis,  nonlinearities,  parameter  uncertainties,  and 
difficulty  in  measuring  variables  such  as  temperature. 

The  control  research  applied  to  SMAs  may  be  divided  into 
three  categories  of  linear  control.  Pulse  Width  Modulation 
(PWM),  and  nonlinear  control.  Different  variations  of  linear 
Proportional  Integral  Derivative  (PID)  controls  have  been 
explored  by  several  researchers  [l]-[6]  while  many  others  have 
used  PWM  [7]— [10].  Several  nonlinear  control  schemes  such 


as  fuzzy  logic,  feedback  linearization,  sliding  mode,  adaptive, 
and  variable  structure  control  have  also  been  explored  by 
researchers  [11  ]— [23] . 

In  this  work,  for  controlling  the  angular  position  of  an 
SMA-actuated  manipulator,  we  have  developed  a  sliding  mode 
controller  that  works  in  conjunction  with  a  classical  feedback 
controller.  The  feedforward  sliding  mode  controller  computes 
the  desired  stress  of  the  SMA  wire  according  to  the  desired 
angular  position  of  the  arm  and  the  feedback  controller  adjusts 
the  control  input  (voltage)  to  maintain  the  desired  stress  in  the 
wire.  The  controller  uses  an  Extended  Kalman  Filter  (EKF) 
to  estimate  the  stress  of  the  SMA  wire.  The  arm  is  actuated 
by  a  bias  type  actuator  constructed  with  SMA  wire,  pulleys, 
and  a  linear  spring.  A  nonlinear  model  is  developed  based  on 
the  arm  nonlinear  dynamic  model,  an  SMA  wire  constitutive 
model,  an  SMA  phase  transformation  model  [24],  [25],  and 
a  nonlinear  heat  convection  model  [14],  [26].  This  model  has 
been  experimentally  verified  by  Elahinia  and  Ashrafiuon  [14]. 
The  Extended  Kalman  Filter  has  been  designed  for  an  SMA 
actuated  manipulator  and  was  tested  through  simulations  in 
our  other  works  [27],  [28].  The  EKF  utilizes  the  control  input 
and  the  angular  position  of  the  arm  to  predict  the  other  state 
variables  of  the  system:  angular  velocity,  SMA  wire  stress, 
and  SMA  wire  temperature.  Since  the  sliding  mode  controller 
transforms  the  nonlinear  dynamics  of  the  arm  into  a  first  order 
system  the  resulting  control  algorithm  has  a  better  overall 
performance  both  in  stabilization  and  tracking  in  the  presence 
of  model  uncertainties. 

II.  The  SMA  Arm  Model 

The  single  degree-of-freedom  (1-dof)  SMA-actuated  arm  is 
shown  in  Fig.  1.  We  have  used  150 jim  diameter  Ni-Ti  SMA 
wire  which  is  actuated  by  electrical  heating  [14].  The  net 
actuating  torque  is  the  difference  between  the  resulting  bias 
spring  and  SMA  wire  torques.  The  SMA  arm  model  consists 
of  phase  transformation,  heat  transfer,  SMA  wire  constitutive, 
and  arm  dynamic  blocks. 

A.  Kinematic/Dynamic  model 

The  nonlinear  dynamic  model  of  the  arm  including  spring 
and  payload  effects  is  represented  by: 

Ie6  =  Tw(cr)  -  rg(0)  -  rs(0)  -  c0  (1) 
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where  9  is  the  angular  position  of  the  arm,  rW)rg,  and  ts 
are  the  resulting  torques  from  SMA  wire,  gravitational  loads, 
and  the  bias  spring,  respectively,  and  cr  is  the  wire  stress.  Ie 
is  the  effective  mass  moment  of  inertia  of  the  arm,  gripper, 
and  the  payload,  and  c  is  the  torsional  damping  coefficient 
approximating  the  net  joint  friction. 

The  SMA  wire  strain  rate  e  and  joint  angular  velocity  9  are 
related  kinematically  as: 


where  rp  is  pulleys  radius  and  Iq  is  the  initial  length  of  the 
SMA  wire. 

B.  Wire  constitutive  model 

The  wire  constitutive  model  shows  the  relationship  between 
stress  rate  (&),  strain  rate,  and  temperature  rate  (T)  [29]: 

a  =  Di  +  0tT  +  fl£  (3) 

where  D  =  Dm+Da  is  (the  average)  Young  modulus,  DA  is 
Austenite  Young  modulus.  Dm  is  Martensite  Young  modulus, 
9t  is  thermal  expansion  factor,  Cl  =  —Dsq  is  phase  transfor¬ 
mation  contribution  factor,  and  e o  is  the  initial  (i.e.  maximum) 
strain  [24], 


Fig.  1 .  The  1  -dof  SMA  arm,  actuated  by  NiTi  wire  and  a  bias  spring 


C.  SMA  Wire  Phase  Transformation  Model 

Due  to  hysteretic  behavior  of  SMA  wire,  the  phase  transfor¬ 
mation  equations  are  different  for  different  load  and  heating 
conditions.  The  reverse  transformation  from  Martensite  to 
Austenite  takes  place  as  a  result  of  heating  or  unloading 

*;  =  -%- sin [aA{T  -  As  ~ -^-)}[aAf  -  ^-a)  (4) 

l  La  Ga 

if  Ca(T  -  Af)  <  a  <  Ca{T  -  As).  Where  0  <  £  <  1 

is  the  Martensite  fraction  coefficient,  £ m  is  the  maximum 

Martensite  fraction  obtained  during  cooling,  T  is  the  SMA 

wire  temperature,  As  and  Af  are  Austenite  phase  start  and 

final  temperatures,  and  aA  =  A  'n_A  ,  bA  =  and  CA  are 

curve  fitting  parameters.  The  forward  transformation  equation 

from  Austenite  to  Martensite  takes  place  because  of  cooling 

or  loading 

£  =  -1  ^  sin [aM{T  -  Mf  -  )}[aMt  -  &\  (5) 

l  G  M  G  M 

if  Cm(T  -  Ms)  <  a  <  Cm{T  -  Mf).  Where  is  the 
minimum  Martensite  fraction  obtained  during  heating,  Ms 
and  Mf  are  Martensite  phase  start  and  final  temperatures, 
and  «M  =  MaZM  •  bM  =  and  Cm  are  curve  fitting 

parameters. 

D.  Heat  transfer  model 

The  SMA  wire  heat  transfer  equation  consists  of  electrical 
heating  and  natural  convection: 

mcpfjf  =  hAc(T  -  T^)  (6) 

where  R  is  resistance  per  unit  length,  cp  is  the  specific  heat, 
to  is  mass  per  unit  length  and  Ac  is  circumferential  area  of  the 
SMA  wire.  Also,  V  is  the  applied  voltage,  is  the  ambient 
temperature,  and  h  is  heat  convection  coefficient. 


Fig.  2.  The  static  stress  of  the  wire  as  a  function  of  the  angular  position  of 
the  arm 


Feedforward  Controller 


Feedback  Controller 


Fig.  3.  The  feedforward-feedback  controller  schematic  and  the  Extended 
Kalman  Filter 
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III.  Control 

The  SMA  phase  transformation  temperatures  are  affected 
by  the  stress  of  the  wire.  Therefore  when  the  arm  moves 
beyond  a  certain  angular  position,  at  which  the  stress  of  the 
wire  is  maximum,  the  transformation  temperatures  decrease. 
Figure  2  shows  the  static  stress  of  the  wire  as  a  function  of 
the  angular  position  of  the  arm.  If  the  desired  angular  position 
is  beyond  this  maximum  stress  position,  the  control  system 
should  first  apply  high  enough  voltage  in  order  for  the  wire 
to  go  through  phase  transformation  and  therefore  reach  the 
maximum  stress.  Beyond  this  position,  while  the  position  error 
is  still  not  zero,  the  voltage  should  drop  rapidly  in  order  to 
prevent  the  wire’s  temperature  from  going  beyond  Af  which 
causes  the  full  phase  transformation.  Therefore,  if  a  control 
system  works  based  on  the  position  error  only,  even  when  the 
arm  is  beyond  the  maximum  stress  position  it  will  apply  high 
voltage  to  the  wire;  hence  the  wire  will  be  overheated  and  the 
arm  will  exceed  the  desired  angular  position.  Although  it  is 
possible  to  tune  such  a  position-based  controller  to  work  well 
for  certain  positions  it  cannot  work  well  for  all  the  set-points. 

One  way  to  solve  this  problem  is  to  design  a  controller 
that  works  based  on  the  stress  of  the  SMA  wire.  In  this 
section  a  sliding  mode  controller  is  designed  to  calculate  the 
desired  stress  of  the  wire  based  on  the  desired  angular  position. 
A  second  feedback  controller,  which  is  a  PID  controller,  is 
used  to  maintain  this  desired  stress  by  applying  the  required 
voltage  to  the  SMA  wire.  Since  this  two-part  controller  tracks 
the  desired  stress  of  the  wire,  when  the  arm  reaches  the 
maximum  stress  position,  the  controller  adjust  the  control 
input  accordingly.  Therefore  the  wire  will  not  overheat.  The 
control  system  is  shown  schematically  in  Fig.  3.  The  actual 
stress  of  the  wire  is  estimated  using  an  Extended  Kalman 
Filter,  which  was  designed  in  our  previous  work  [27]. 

Sliding  control  is  a  robust  nonlinear  Lyapunov-based  control 
algorithm  in  which  an  nth  order  nonlinear  and  uncertain 
system  is  transformed  to  a  1st  order  system.  Let’s  consider 
the  SISO  system 

x(n)  _  _|_  5(x)m  (7) 

where  u  and  x  are  the  scalar  input  and  output  respectively  and 
x  is  the  state  vector.  For  this  system  the  nonlinear  functions 
/  and  b  are  not  exactly  known.  The  control  objective  is  for 
the  system  to  follow  a  desired  state  vector  x^.  Furthermore 
consider  the  time-varying  surface  S(t)  in  the  state-space 

s(x,f)  =  (^  +  A)n_1S  (8) 

where  x  =  x  —  Xd  is  the  tracking  error  and  A  is  a  positive 
constant.  It  can  be  shown  that  tracking  the  desired  state  vector 
is  equivalent  to  remaining  on  the  surface  S(t)  [30].  Therefore 
the  control  law  should  be  selected  in  a  way  that  the  distance 
to  this  surface  decreases  along  all  system  trajectories 


where  77  is  a  positive  constant;  equation  9  is  called  the  sliding 
condition.  As  the  result  all  the  trajectories  are  forced  to  reach 
the  sliding  surface  in  finite  time  and  stay  on  the  surface  for  all 


future  times.  S(t)  is  called  the  sliding  surface  and  the  systems 
behavior  once  on  the  surface  is  called  sliding  mode  (s  =  0). 
It  can  be  shown  that  the  reach  time  from  the  initial  state  to 
the  surface  is  [30] 

breach  <  s(t  =  0)/t]  (10) 

Thus  a  typical  motion  under  sliding  mode  control  consists  of 
reaching  phase  and  sliding  phase  during  which  the  motion  is 
confined  to  the  sliding  surface.  Further  details  on  the  sliding 
mode  control  can  be  found  in  [30],  [31]. 

A.  Controller  design 

The  sliding  controller  in  this  work  is  deigned  to  calculate 
the  desired  stress  of  the  SMA  wire.  Considering  the  equation 
of  motion  (equation  1),  the  sliding  surface  can  be  defined  as 

S=(ft+X)9  (11) 

where  9  =  9  —  9d-  The  dynamic  of  the  system  once  on  the 
sliding  surface  can  be  written  as 

s  =  0  (12) 

This  equation  describes  the  dynamic  of  the  system  while  in 
the  sliding  mode  and  can  be  used  to  find  the  desired  torque 
of  the  wire: 

s  =  9  -  9d  +  X9  =  (rw  -  rg  -  rs  -  c9)/Ie 

-  6d  +  \e  =  0  (13) 

and  therefore 

Tw  —  Tg  +  Ts  +  c9  T  Ie9d  IeX9 

+  Afsat(s)  (14) 

where  K  is  the  control  gain.  The  saturation  term  is  added 
to  the  desired  torque  in  order  to  address  the  uncertainties  in 
the  model.  This  term  also  causes  the  trajectories  to  merge  to 
the  sliding  surface  when  the  initial  state  of  the  system  is  not 
on  the  sliding  surface  (0(0)  ^  0).  The  desired  stress  of  the 
SMA  wire  can  be  readily  calculated  using  the  desired  torque 
of  equation  13  . 

IV.  Results 

The  nonlinear  SMA  model  used  for  the  simulations  in  this 
work  has  been  previously  verified  through  several  experiments 
[14],  Figure  4  shows  the  simulation  and  experimental  results 
for  two  different  voltages.  Using  this  model,  in  this  section,  the 
simulation  results  with  the  sliding  mode  controller  (SLD)  is 
presented  and  compared  with  a  PID  controller.  Figure  5  com¬ 
pares  the  regulation  performance  of  the  proposed  controller 
with  that  of  a  PID  controller.  The  PID  controller,  which  works 
based  on  the  position  error,  can  be  tuned  to  perform  well  for 
certain  set-points  while  the  same  SLD  controller  works  well 
for  all  the  desired  positions.  The  regulation  performance  of 
the  SLD  controller  for  several  set-points  is  shown  in  Fig.  6. 
It  can  be  seen  that  the  arm  reaches  the  desired  position  for  all 
the  set-points  in  less  than  4  sec  even  though  small  overshoots 
exist.  The  performance  of  the  controller  in  tracking  the  desired 
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Fig.  4.  Comparing  open-loop  simulation  and  experimental  results 


Time(sec) 

Fig.  7.  The  tracking  performance  of  the  feedback  part  of  the  controller 
following  the  desired  stress  of  the  wire  6 d  =  45° 


time(sec) 


Time(sec) 


Fig.  5.  Comparing  the  regulation  performance  of  the  feedforward-feedback  o  T,  ,  ,•  r  f  ..  f  ,r  ,  r  ,,  ,  .  ,, 

I.  •,  rig.  o-  Ihe  tracking  performance  or  the  reedrorward-ieedback  controller 

controller  (SLD)  with  a  PID  controller  (SLD)  following  a  changing  step  input 


Fig.  6.  The  regulation  performance  of  the  feedforward-feedback  controller  Fig  9  xhe  pacing  performance  of  the  feedforward-feedback  controller 
(SLD)  for  desired  positions  0^  =  —40°  up  to  0^  =  60°  (SLD)  following  a  repeating  step  input  changing  between  0°  and  45° 
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stress  is  shown  in  the  Fig.  7.  The  tracking  performance  is  good 
and  as  the  result  arm  reaches  the  desired  angular  position. 

Figures  8,  9,  and  10  show  the  tracking  performance  of 
the  controller  for  repeating  step  and  ramp  inputs  respectively. 
The  arm  closely  follows  the  desired  angular  position;  however 
because  of  the  speed  limitation  of  the  SMA  wire,  the  tracking 
performance  is  better  when  arm  follows  desired  position 
signals  with  smaller  motion  spans.  Figure  1 1  and  12  show  the 
tracking  performance  of  the  controller  for  sinusoidal  inputs 
with  two  different  frequencies.  Although  the  motion  span  is 
the  same  for  the  two  sinusoidal  desired  inputs,  it  can  be  seen 
that  the  tracking  performance  is  better  for  slow  varying  inputs 
due  to  inherent  speed  limitation  of  the  SMA  wire. 

Figure  13  shows  the  effect  of  the  sliding  surface  parameter 
A  on  the  regulation  performance  of  the  controller.  The  sliding 
mode  control  transforms  the  dynamic  of  the  system  into  a  1st 
order  system  s  =  ( +  A )9.  Therefore  increasing  A  is  similar 
to  adding  more  damping  to  the  system  hence  the  overshoot 
decreases  and  the  system’s  response  gets  slower.  Figure  14 
shows  the  effect  of  discontinuity  gain  K  on  the  regulation 
performance  of  the  controller.  Increasing  this  gain  is  equivalent 
to  putting  more  energy  to  the  system  and  therefore  the  system 
responds  faster  and  with  higher  overshoot;  decreasing  the  gain 
has  the  opposite  effect.  Thus  the  performance  of  the  controller 
can  be  augmented  by  adjusting  the  gain  K  and  weighting 
factor  A  of  the  sliding  mode  controller. 

Figures  15  and  16  depict  the  robustness  of  the  controller. 
The  effective  moment  of  inertia  and  the  damping  coefficient 
are  two  of  the  parameters  that  are  used  both  in  the  dynamic 
model  and  the  sliding  mode  control  law  (equations  1  and 
14).  To  model  uncertainties  different  values  are  used  for  these 
parameters  in  the  sliding  mode  controller  and  in  the  model. 
From  equation  14,  it  can  be  seen  that  if  the  moment  of  inertia 
used  by  controller  Ie  is  larger  than  the  actual  moment  of 
inertia,  system’s  response  is  more  damped.  On  the  other  hand 
if  the  damping  coefficient  used  by  the  sliding  mode  controller 
c  is  larger  than  the  actual  c,  the  system  responds  similar  to 
a  less  damped  system  and  therefore  with  larger  overshoot. 
Increasing  damping  coefficient  results  in  calculating  a  larger 
desired  stress,  the  feedback  controller  applies  larger  inputs 
to  the  wire  and  therefore  there  exists  larger  overshoot  and 
faster  response.  In  the  same  way,  choosing  a  larger  moment 
of  inertia  reduces  the  calculated  desired  stress  and  therefore 
the  controller  applies  smaller  voltage  to  the  wire  and  hence 
the  system  responds  slower  and  with  smaller  overshoot. 

V.  Conclusion 

A  sliding  mode  controller  is  designed  for  a  rotary  SMA- 
actuated  manipulator  and  has  been  tested  through  simulations. 
This  control  law  solves  one  of  the  challenges  for  control  of  this 
system  (i.e.  dependency  of  the  transformation  temperatures 
on  the  wire  stress).  The  controller  regulates  the  wire’s  input 
voltage  based  on  the  desired  stress  calculated  by  the  sliding 
mode  controller.  An  Extended  Kalman  Filter  is  used  for 
estimating  the  SMA  wire’s  stress.  Satisfactory  control  results 
have  been  shown  in  terms  of  the  reaching  desired  angular 
position  and  tracking  desired  trajectories.  Since  the  model 


Time(sec) 

Fig.  10.  The  tracking  performance  of  the  feedforward-feedback  controller 
(SLD)  following  a  repeating  step  input  changing  between  —45°  and  60° 


Time(sec) 

Fig.  11.  The  tracking  performance  of  the  feedforward-feedback  controller 
(SLD)  following  a  sinusoidal  input  with  frequency  uj  =  0.6rad/sec 


Fig.  12.  The  tracking  performance  of  the  feedforward-feedback  controller 
(SLD)  following  a  sinusoidal  input  with  frequency  u ;  =  2.0rad/sec 
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Fig.  13.  Effect  of  the  sliding  surface  parameter  A  of  the  regulation 
performance  on  the  controller,  K  =  1 


Fig.  14.  Effect  of  the  discontinuity  parameter  across  the  sliding  surface  K 
on  the  regulation  performance  of  the  controller,  A  =  30 

had  been  previously  verified  against  experimental  data  and 
the  controller  is  designed  based  on  this  model,  it  is  expected 
that  the  controller  works  well  with  the  SMA  actuated  arm.  It 
is  also  shown  that  the  control  performance  is  robust  in  the 
presence  of  model  uncertainties. 
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